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Results for the production of W^-bosons in conjunction with up to three jets including parton 
shower corrections are presented and compared to recent LHC data. These results consistently 
incorporate the full next-to leading order QCD corrections through the Mc@NLO method, as imple- 
mented in the SHERPA event generator, with the virtual corrections obtained from the BLACKHAT 
library. 



The production of VF-bosons accompanied with QCD 
jets is a Standard-Model reaction central to the exper- 
imental program at the LHC, as it constitutes an im- 
portant background to many new physics searches. It is 
also of interest in its own right, to study jet production 
and evolution in a hadron-collider environment [TJ [2] , to 
improve the jet energy scale determination by the exper- 
iments, or to study multiple parton scattering processes. 
At the LHC as well as at the Tevatron, typically, good 
agreement is found when comparing respective data with 
next-to-leading order (NLO) perturbative QCD predic- 
tions, which have recently been made for PF-boson pro- 
duction in association with three and even four jets [3HS]. 

The development of two methods to match NLO ma- 
trix elements with parton showers, MC@NLO [6J and 
POWHEG [7], allowed, at the same time, to include earlier 
parton-level calculations into parton-shower Monte-Carlo 
(PSMC) programs, among others this includes Z+l-jet 
production in the POWHEGBOX |8J and VF+2-jets pro- 
duction in aMc@NLO [9]. In this letter we present the 
implementation of the VK+jets process in a variant of the 
MC@NLO method [TO], including up to three jets at NLO 
accuracy. Our implementation in SHERPA [TJJ can eas- 
ily be extended to other processes with similar or even 
higher final state multiplicities. For the process at hand, 
we use the BLACKHAT library for NLO calculations to 
evaluate the finite part of virtual corrections [3]. 

Our new scheme to implement the Mc@NLO technique 
is based on the exact exponentiation of Catani-Seymour 
dipole subtraction terms [TO]. This method allows to 
circumvent the otherwise occurring integral over resid- 
ual real-radiative contributions to the NLO cross sec- 
tion, that arise from the modified subtraction scheme in 
MC@NLO [5J. It also allows, for the first time, to obtain 
the correct soft-gluon limit in the first emission of the 
PSMC, such that no adjustments to the splitting kernels 
of the parton shower must be made in the soft region 
to match the correct behavior of the subtraction terms. 
Due to the existence of multiple color structures at Born 



level, this point is particularly important when simu- 
lating processes with Mc@NLO techniques that involve 
a large number of QCD partons, like W+3 jets. Our 
method allows to prove that the approximations which 
are usually made in Mc@NLO to deal with soft diver- 
gences are justified. In fact it can be shown to correctly 
take into account the full color structure of the processes 
at NLO. 

The Mc@NLO cross section can be written in the 
form [5JCEU] 
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The terms B, V, V- s \ and D^ s ^ represent the Born con- 
tribution, virtual correction plus collinear counterterms, 
integrated subtraction terms and real subtraction terms. 
3>b and denote Born- and real-emission phase space 
with = <E , B<g>$i. Real-emission matrix elements R are 
separated into an infrared-singular (soft) and an infrared- 
regular (hard) part, D^ A ) and H< A ) as R = D< A ) + H^ A ). 
This leads to the definition of the Sudakov form factor 



A( A >(M') = 
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The key point of our new technique is that the integral 
in Eq. (2]) is avoided if D (A) = D (s) , i.e. if the subtrac- 
tion kernels are employed for parton showering. This 
can be achieved using Catani-Seymour subtraction, by 
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dynamically correcting a parton shower based on spin- 
and color-averaged splitting operators. The method was 
applied previously to the W ± - Z- and Higgs+l-jet pro- 
duction processes [TO]. In this publication we show that 
it is not limited to the case of one final-state parton at 
Born level. We present results for W+2- and W+3-jet 
production, which contain the most general color struc- 
tures. The latter of these two processes has neither been 
implemented in the Mc@NLO nor in the POWHEG tech- 
nique so far. 

We use the framework of the SHERPA event genera- 
tor [TTJ, including its automated Mc@NLO implementa- 
tion [TU], the matrix-element generator AMEGIC++ [T2"] . 
an automated implementation |13| of the Catani- 
Seymour dipole subtraction method [14] and the parton 
shower model described in |15j . This parton shower uses 
transverse momentum ordering. We thus avoid the prob- 
lem of truncated emissions [7J. We restrict the resum- 
mation region as described in [TD] by setting a = 0.01. 
We appreciate that this is a process-dependent param- 
eter, which in the future is going to be replaced by a 
phase-space cut in terms of transverse momentum. Such 
a choice is physically more sound and will ultimately 
allow for a fully automated implementation which cor- 
rectly resums all leading logarithms. The finite part of 
virtual corrections is computed using the BLACKHAT li- 
brary [4j. For the W+3-jet virtual matrix element we use 
the leading-color approximation in BLACKHAT to avoid 
an unnecessary increase in CPU time for the simulation. 
Subleading color configurations in virtual corrections of- 
ten play a minor role in W+multi-jet processes |16j . 
They might, however, be important in other situations. 
As we focus on the interface between the NLO calcu- 
lation and the parton shower in this publication these 
effects are neglected. The CTEQ6.6 PDF set [T7j is em- 
ployed together with the corresponding parametrization 
of the running coupling. Following |18j renormalization 
and factorization scales are chosen as fj,R = fj,p = H' T , 

where H' T = l/2^Jj2PT,j + E'tw Predictions are pre- 
sented at three different levels of event simulation: 

"NLO": Fixed-order matrix-element calculation, 

"MC@NLO lem": MC@NLO including only the first 
(hardest) emission in the parton shower, 

"MC@NLO PL": Mc@NLO including full parton 
showering, but no non-perturbative effects. 

Non-perturbative effects are neglected in this studjF] in 
order to focus on the features of the Mc@NLO method in 
processes with intricate color topologies. The effects of 



1 The observables displayed here are relatively insensitive to non- 
perturbative corrections. 



multiple parton interactions, hadronization, hadron de- 
cays and final state QED radiation have been analyzed 
in detail in [TDJ. Likewise we do not quote scale uncer- 
tainties of the NLO results. Those have been analyzed 
in great detail in [5] and [TU]. The aim of this study is 
instead to present an application of the Mc@NLO variant 
suggested in [TD] to complex QCD final states in order to 
verify its universality and test its versatility. 

The analysis is carried out with the help of Rivet [Jj5] 
following a recent study of W r± +jets production by the 
ATLAS collaboration [2]. Events are selected to contain 
a lepton within \rj\ < 2.5 withp^ > 20 GeV and requiring 
E™ iss > 25 GeV. A cut on > 40 GeV is additionally 
applied. All particles other than the leading electron 
and neutrino are clustered into anti-fct jets with R = 
0.4 and p± > 30 GeV. The analysis is carried out in 
jet multiplicity bins up to N = 3 and cross sections are 
studied differentially in several observables. 

The results for each observable are predicted at NLO 
accuracy, i. e. all differential cross sections for Vl /± + > n- 
jet events are generated using the W ± + n-jet NLO or 
MC@NLO calculation. For n > 0, the VF+n-jet matrix 
element must be regularized by requiring at least n jets 
with a minimum transverse momentum. This cut is cho- 
sen to be p^ n > 10 GeV to make the event sample inclu- 
sive enough for the analysis. We have checked that our 
results are independent of the precise value of this cut by 
varying it from 5 to 15 GeV in every individual jet bin. 

Table [T] compares total cross sections in four inclusive 
jet multiplicity bins. The ATLAS measurement is repro- 
duced very well both by the fixed order calculation as well 
as by the Mc@NLO matched simulation. The agreement 
between the NLO results and the MC@NLO simulation 
is excellent, indicating that the matching to the parton 
shower does not alter the jet production rate as predicted 
by the fixed-order calculation. 

In Fig. [I] we display a comparison of the transverse 
momentum spectra of the first, second and third hardest 
jet in W+ > 1-, 2- and 3-jet production. No significant 
changes are observed when switching from the fixed-order 
calculation to the Mc@NLO simulation, again indicating 
that the hard kinematics predicted by the NLO result are 
respected in the subsequent parton-shower evolution. 

Fig. [2] focuses on W+ > 2-jet events. Angular correla- 
tions between the two leading jets are sensitive to QCD 
corrections in the W + 2-jet process and are thus a useful 
observable to validate the QCD radiation pattern which 
is generated in our Mc@NLO. Both, the rapidity and 
azimuthal separation of the jets are predicted in perfect 
agreement with data. 

In summary, we have shown in this letter how our re- 
cently proposed method for implementing Mc@NLO with 
the help of Catani-Seymour subtraction and a parton 
shower based on Catani-Seymour subtraction kernels can 
be used to produce novel and relevant results for one of 
the most challenging collider signatures to date. Pre- 
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W ± + > n jets 


ATLAS 


NLO 


Mc@NLO lem 


McQNLO PL 


n = 


5.2 ±0.2 


5.06(1) 


5.09(3) 


5.06(3) 


~* i ~* ^ on /" ■ t 

n = 1, > 20 ijev 
p ±J > 30 GeV 


0.9o ± 0.10 
0.54±0.05 


o.yoo(oj 

0.527(4) 


0.968(10] 

0.534(7) 


0.889(10] 
0.474(7) 


Tl — Z, pj_ j ,> ZU Vjre V 

p ±J - > 30 GeV 


U.ZO zt U.U4t 

0.12±0.02 


U.ZOO^Z ) 

0.120(1) 


U.ZDU^O ) 

0.123(2) 


u.zooi ^ ^ 
0.109(2) 


rc = 3, p x i > 20 GeV 
p ±j > 30 GeV 


0.068 ±0.014 
0.026 ±0.005 


0.072(3) 
0.026(1) 


0.059(3) 
0.022(2) 


0.060(3) 
0.021(1) 



TABLE I. Total cross sections in nb for W ± ± > 0, 1, 2, 3 jet production as measured by ATLAS [2] compared to predictions from 
the corresponding fixed order calculations, and matrix-element/shower level Mc@NLO simulations. Statistical uncertainties of 
the theoretical predictions are quoted in parentheses. 



dictions for W±3-jets production were made, a process 
which has not been attempted in either the POWHEG 
or the MC@NLO framework so far. We have compared 
results for W+0-, 1-, 2- and 3-jet production to recent 
ATLAS data and found excellent agreement for all ob- 
servables, with only a selection of them presented here. 
The success and the simplicity of our Mc@NLO vari- 
ant make it a prime candidate for the implementation 
of a matrix-element parton-shower merging algorithm at 
next-to-leading order. 
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FIG. 1. Transverse momentum of the first, second and third 
jet (from top to bottom) in W + > 1, 2, 3 jet production as 
measured by ATLAS [2] compared to predictions from the 
corresponding fixed order and Mc@NLO simulations. The 
gray bands display the combined statistical and systematic 
experimental uncertainties. 
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FIG. 2. Angular correlations of the two leading jets in 
W + > 2 jet production as measured by ATLAS [2] com- 
pared to predictions from the W + 2 jet fixed order and 
Mc@NLO simulations. The gray bands display the combined 
statistical and systematic experimental uncertainties. 



